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Abstract 

Starting with a Chern-Simons theory, we derive an effective action for inter- 
acting quantum Hall skyrmions that takes into account both large-distance 
physics and short-distance details as well. We numerically calculate the classi- 
cal static skyrmion profile from this action and find excellent agreement with 
other, microscopic calculations over a wide range of skyrmion sizes including 
the experimentally relevant one. This implies that the essential physics of 
this regime might be captured by a continuum classical model rather than 
resorting to more microscopic approaches. We also show that the skyrmion 
energy closely follows the formula suggested earlier by Sondhi et al. for a 
broad parameter range of interest as well. 
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I. INTRODUCTION 

Quantum Hall systems! near filling factor v = l/(2n + 1) should manifest a charged 
topologically stable object called a skyrmion. Various experiments have shown strong 
evidence for the existence of this exotic object in two dimensional electron systems. In 
particular the ground state of the quantum Hall system near filling factor z/ = 1 is believed 
to be a manybody state of weakly interacting skyrmions.&0 

Just as the thermodynamics of certain superconductors can be well-described by their 
vortex degrees of freedom alone, skyrmions, the defects of the incompressible quantum Hall 
liquid, may well describe the essential physics of the quantum Hall system near u = 1. The 
thermodynamics of vortex systems can be derived from the phenomenological Ginzburg- 
Landau (GL) description of superconductivity. In traditional BCS-theory the spin-singlet 
Cooper-pairing mechanism removes the electronic spin degrees of freedom from the dynam- 
ics. Hence vortices have no spin structure,lll and are entirely described by their integer 
charge. In contrast, the quantum Hall skyrmion is a topological spin-texture made of an 
intricate pattern of electron-spin orientations. It possesses internal degrees of freedom asso- 
ciated with this texture that allow for novel dynamics.0 

In analogy to GL theory in the superconductivity, the Chern-Simons theory of quantum 
Hall systems was derived by Zhang et a/.0 and subsequently extended by Lee and Kandl to 
describe the possible spin-unpolarized quantum Hall liquid. In superconductivity, it has been 
very useful to obtain a phenomenological theory of vortices based on GL theory.0 A similar 
theory for skyrmions should be invaluable for understanding transport, phase transitions, 
and spatial ordering.^ 

In this paper, starting from the Chern-Simons theoryi'0 we derive an action for the 
many skyrmion system which takes into account both large-distance physic^ and relatively 
short-distance details as well. Using this action we perfrom numerical studies of the classical 
skyrmion solution. Recently Abolfath et al^ pointed out for small skyrmions, i. e. those 
typical in GaAs samples,0'S0 the minimal effective field theory does not give a good quan- 
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titative agreement with the Hartree-Fock or exact- diagonahzat ion study. We find that the 
various properties of a static skyrmion solution obtained from the above action exhibit an 
excellent agreement with the microscopic study for 'baby' skyrmions (textures containing 
from two to about 10 fiipped spins) and larger. This implies that current experiments can be 
well modelled by a continuum theory, obviating the need for the more involved microscopic 
approaches. 

In section II, we briefiy summarize the Chern-Simons description of quantum Hall effect 
and then derive an action for many skyrmion system. In section III, the classical skyrmion 
profile is numerically calculated. The various properties of a static skyrmion solution are 
compared with the microscopic study. We show that the skyrmion energy follows the formula 
suggested earlier by Sondhi et al. very accurately for a broad parameter range. In section 
IV, we conclude with a summary. 



II. CHERN-SIMONS THEORY 

We briefiy summarize the Chern-Simons description of quantum Hall effect which was 
subsequently extended by Lee and Kane in order to incorporate a possible spin-unpolarized 
quantum Hall liquid.01i In the bosonic Chern-Simons theory0, an electron is viewed as 
a composite object of a boson and a fiux tube carrying an odd-multiple of fiux quanta 
(po = h/e attached via a Chern-Simons term, which correctly ensures fermionic statistics for 
the electron. We begin with the effective Chern-Simons Lagrangian introduced by Lee and 
Kanei 



a^] = *^(9o - mo)^a + ^ 



,V 
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+ \jdr' ('l'.(r)^.(r) - p) V{\r - r'|) ('l'.(r')^.(r') - p) + C,, (2.1) 

where represents a bosonic field with spin cr = ±i, m the effective mass of boson, 
Aex a vector potential for the external magnetic field, p a mean density of boson, and 
V{r) = e^/er Coulomb interaction between electrons. We will use the convention of units 
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where h = c = e = 1. The Chern-Simons term Ccs can be written as 

-^cs = -^;^^fiuxci'fj.duax (2.2) 

where is a statistical gauge field with yU = 0, 1, 2, (1 and 2 are spatial indices, is time) 
and we impose a Coulomb gauge V ■ a = 0. The variable a is taken to be an odd integer, 
a = 2n + 1, in order to describe a fermionic system. 

In order to separate the amplitude and the spin degrees of freedom from the bosonic field 
we introduce the CP^-field Za-, which satisfies \z„\ = 1 and is related to by 



+ 2 



K V - 



Z(j . Zcr a. Ag 



= ^Jo z^{x) (2.3) 

where Jo represents the boson density. The Lagrangian in terms of the fields Jo and z^ is 
given by 

1 r d 
C[z„, a^] = - J dr' (Jo(r) - p) V{\v - r'|) (Jo(r') - p) + iJo {z^-t-z„ - ao) 

C^La + C,, (2.4) 

where k = p/m and the gradient term in Jo(r) is neglected.0 The local spin orientation m 
is related to the field z^^ via m = ZaaapZp, so that the static non-linear sigma model (NLa) 
term can be written by 

£nl. = ^(Vm)2. (2.5) 

o 

In order to decouple the quartic term, one introduces a Hubbard-Stratonovich field which 
represents a bosonic current J(x). After integrating out the statistical gauge field a^, one 
obtains the following LagrangianiUll 

C[z^, A,] = \j dr' ( Jo(r) - p) V{\r - r'|) ( Jo(r') - p) + i-| Jp + C^^^ 

+ t{A, + 4°)) j; - '-aA, (J, - JW) + 1-Ao {B - 2'nap) (2.6) 

where J^ = J^^^ + ■^t^j.xdvAx and J^°^ = (p, 0, 0) = ^t^y\dyA^^\ Here we used the fact 
that the bosonic current and density satisfies a continuity equation 9^ J^ = 0. The skyrmion 
three-vector Jf, can be written in terms of m: 



Ja = ^t^^x^dyva. X ^Ain) ■ m (2.7) 

where we explicitly put the factor of 1/a in the definition of J*. The Jq{x) is the topological 
charge density of the spin texture, which is proportional to the electronic charge density. 

Zhang et al. have shown that when external magnetic field is tuned so that the number 
of flux quanta is commensurate with the mean boson density, {B = 27cap, so that, u = 
l/(2n + 1)), the bosons condense and form a superfiuidlli. The bosonic superfluidity in the 
Chern-Simons theory implies that a quantum Hall effect occurs in the corresponding two- 
dimensional electron system. The ground state has been shown to be a fully spin-polarized 
quantum ferromagnet. By noticing that the dynamics of quantum Hall system with a spin- 
polarized ground state will follow that of a quantum ferromagnet and that the skyrmion is 
a charged object of the system^, Sondhi et al. proposed a phenomenological action, which 
is valid for the long- wavelength and small-frequency limit. 

However we can explicitly integrate out the bosonic field in Eq. ( p.6| ) and derive an action 
for skyrmions which takes fully into account the short- and long-distance physics. A similar 
exercise is standard in the GL theory for vortices.0ii We proceed by integrating out the 
bosonic field with a Coulomb gauge condition, V- A = 0. In order to impose the Coulomb 
gauge, we introduce an auxilliary field A(x) and introduce an additional term iX{x)V -A into 
Eq.( p.6|) . Since the action is still quadratic in the field and diagonal in the frequency- 
momentum space, one can exactly integrate out the bosonic field A^. The integration of 
the k = mode leads to the following relations between the skyrmion density Ug and the 
external magnetic field B, i.e. Ug = {27Tap — B)/ (27r). 

After lengthy but straightforward calculations, we finally obtain the following action 
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SE[ni] = ^Y.Vik) |J„^|2 + ^1 |P|2 + IdTjdr £nl<x + J dr J dr 



where k stands for (cj,k), P{k,uj) = na'^ + V{k)k'^/{A7r'^) +uj'^/{4:7r'^K), and = acj/(27r) 
and the Zeeman term is given by 

£. = ^ (l-m.(r)) (2.9) 

where the magnetic length i = {hc/\e\BY^'^ and t = (l/2)giXBB. In the above action, 
the first term represents a charge-density interaction between skyrmions including the self- 
energy contribution. The function V{k) is the Fourier transform of the Coulomb interaction; 
due to the additional term V{k)k'^/ (47r^) oc A; in the denominator, the interaction is modified 
by the short-range fluctuations of the gauge field from a Coulombic one ~ 1/r to ln(l/r) at 
short-distances. The next term represents the kinetic energy for skyrmion. In the limit of 
long-wavelengths and small-frequencies, we obtain the following action: 

^EnA:)|Jol^ + E ^ 
+ taJ2A'^''\-k)-3'{k)-J2'^ J^{-k)z-kx 3%k). (2.10) 

The term iaA.^^\vi) ■ J| indicates that the skyrmion views the original boson as a magnetic 
flux tube§. The last term in Eg. (|2.10|) contains the exchange-statistics of skyrmion. It can 
be re-written into the more suggestive form zA^fc ■ J* where 

V X A,fc = 27raJo'. (2.11) 

In order to see the exchange-statistics of a skyrmion, suppose that all the other skyrmions 
are at rest while one moves around the static skyrmion configuration. When a skyrmion 
traverses around a closed loop, this term generates a phase proportional to the number of 
skyrmions enclosed in the loop: 

j dr dt A,k ■ J" = 27ragsk J dr (2. 12) 

where S stands for the space enclosed by the closed skyrmion loop. Using the fact that the 
skyrmion charge gsk is equal to e/a, one can show that skyrmion picks up a phase {27i/a)Ne^c, 
where iVenc is the number of skyrmions enclosed. Now consider a process which exchanges 



two skyrmions: in the rest frame of one of the skyrmions, the exchange corresponds to the 
other skyrmion moving about the first in a half circle; and hence it picks up a phase vr/o;. 
Since a is (2n + 1), the statistical phase of a skyrmion is Tij{2n + 1). For n = 0, skyrmion 
is a fermion, while for 7^ 0, it's an anyon^il. 

III. NUMERICAL SOLUTION OF CLASSICAL SKYRMION 

We calculate the classical skyrmion solution with varying Zeeman energy and make a 
comparison to the microscopic result obtained by Hartree-Fock and exact diagonalization 
studies. We begin with the static energy functional E\m\ derived from Eq.( p^.8| ) 

Note that the charge-density interaction changes from a Coulombic one ~ 1/r to ln(l/r) 
at short-distances. Since the skyrmion size is determined by balancing the Zeeman energy 
and the Coulomb interaction, the size of skyrmion will be reduced from the estimates of the 
minimal field theoryEl. 

By using the fact that the skyrmion solution is azimuthally-symmetric, we choose the 
form m = (sin^(r) cos0,sin^(r) sin0, cos^(r)). In order to solve for the classical skyrmion 
profile, we extremize the energy functional E\m\ with respect to Qir) and obtain the following 
equationlli 

\d sin2g /^^gi^g I ^ /^sing^^^^_^ 

r dr \ dr J 2r^ V n \ n r 

where g = 2t/(e^/e£) and /(r) is given by 



fir) 



I dr' (^-^^U{r,r')yme{r')^ (3.3) 



where the azimuthally-averaged interaction potential, U{r,r'), is given by U{r,r') 



4:Jo^^d(P \/(^(r-r')2 + 4rr'sin2 0) with V{r) = dq Mqr)/{1 + 6g), h = 2^2/^, and 
Jo the Bessel function of the zeroth order. If instead we set 6 = 0, we recover the case of a 
pure Coulomb interaction, and then Eq. (|3.2| ) agrees with the one obtained by Abolfath et 



alx^. We have used the fact that k is equal to 4ps where ps = e^/(16-\/27re£) as shown by 
Sondhi et a/.i and by Moon et a/.i. We impose the boundary conditions: 6{r = 0) = n and 
e{r -4 oo) = 0. 

A brief explanation of how the non-local term /(r) is handled is appropriate. In order to 
calculate the function /(r), we first need to obtain the explicit form of a modified skyrmion 
interaction V{r) following the integration over momenta q by standard numerical integration 
methods. The interaction potential V{r) varies as 1/r for r ^ b and (1/6) ln(l + b/r) 
for r <^ b. By virtue of this asymptotic behaviour, we accurately approximate V{r) by 
(1/6) ln(l+6/r)+l/(ao+air+a2r^-|-a3r^), where the best fit parameters are ao = 13.222, ai = 
6.158,02 = 1.223, and = 0.0004. After differentiating V {^J {r - r'f + Arr' s\n^ 0) with 
respect to r, the function R{r,r') = dU{r,r')/dr can be written as follows 

/■'^/^ , r — r' + 2r' sin^ (t> , , , . 

R{r, r' = 4 / d(^-. , ^ , . 2 ^ H{x) 3.4 

^ Jo (r - r')2 + 4rr' sm^ 

where x = \J [r — r'Y + Arr' sin^ and //(x) = x{dV/dx) satifies lim^_>o -^^(2;) = —1/6. 
We can now perform the integration over 0. By noticing that for r = r', then i?(r, r') as 
a function of r' asymptotically approaches a step-function — (7r/4r)sgn(r' — r), we se that 
i?(r, r') can be decomposed into a regular and a discontinuous part: 

i?(r, r') = ^e(r - r') + /?reg(r, r') (3.5) 

where 6(x) is the Heaviside step function and the R^^g is a smooth, continuous function. 
Then /(r) can be decomposed as well yielding 

7r d 
f(r) = —n + cos e(r)) + dr'R,eJr,r')— cose (r). (3.6) 
2r JO dr 

The integration in the second term can be accurately done by standard numerical methods. 

In order to solve the differential equation numerically, we discretized the equation by Nq 
segments with uniform spacing Ar in units of the magnetic length i. First, the differential 
equation is solved for a typical experimental value of ^ = 0.015, which gives a suitable set of 
parameters for the number of discretization Nq and the spacing Ar. The set of parameters 
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are chosen to be A^o = 500 and Ar = 0.1; the results are not sensitive to reasonable choice 
of Nq and Ar. For other values of g, we re-scale Ar by the length scale ^ oc 1/v^ set by the 



Zeeman energy, which controls the size of skyrmion : Ar = O.lyO.015/^. The differential 
equation is solved using an iterative method on the finite number of grid Nq. 

We first calculate the number of spin-flips as a function of g which is defined as foUowJli 




In Fig.(|lD, the sohd curve is our field-theoretical result for K > 2, which shows a very good 
agreement with the Hartre-Fock calculation by Fertig et a/.0. Since the Hartree-Fock and 
exact diagonalization method can not easily calculate the single skyrmion energy due to 
the long-range nature of Coulomb interaction, the energy cost for creating a charge-neutral 
object is calculated by Abolfath et al^. However it has been pointed out by Abolfath et 
al. that the energy difference A(ii') between the skyrmion with K and K + 1 flipped spins 
can be still obtained. ^{K) corresponds to 2t at which an energy level crossing occurs 
between a skyrmion with K spin-flips and K + 1 spin-flips.El Since K is not quantized in 
the continuum field theory, ^{K) roughly corresponds to the value of 2t at the field where 
the number of flipped spins equals K + 1/2. In Fig.(^, A(_ft') is plotted with respect to K. 
The solid curves are obtained from our field-theoretical calculation. Our field-theoretical 
calculation gives good overall agreement with the microscopic calculations ranging from 
'baby' skyrmions {2 < K < 10) to full-fledged skjnrmions. For 'infant' skyrmions, i.e. K < 2, 
the quantum fluctuations about the classical skyrmion solution will become importantlli. 
Since the number of spin flips in typical experiments with GaAs sample£ll0@ is ~ 3, our 
field-theoretical calculation produces a reasonable result over the whole parameter range of 
interest. For relatively large ('adult') skyrmions, Sondhi et al. have suggested the following 
asymptotic formula for the skyrmion energy B : E[g]/{ATip,) = 1 + A/(47r)[^ln(^)]V3 with 
A = (37r^/4)(72/7r)^/^ = 24.9. One can see that in the absence of Zeeman energy, i.e. 
g ^ 0, the skyrmion energy correctly approaches to the NLcr model result. The above 
formula is not completely on rigorous footing yet as commented by Sondhi et al. in their 






(3.7) 
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recent publicationEll. In Fig.(|]), the skyrmion energy is plotted with respect to [^ln(^)]^/^, 
which shows a nice hnear behaviour for small g. We estimate the line slope by a linear fit 
and A is obtained to be 24.9 ± 0.1. Hence we have numerically confirmed that the skyrmion 
energy follows the formula quite accurately over a very broad parameter range. For example, 
for typical experimental value of ^ = 0.015 where K ~ 3.7, the energy obtained from the 
above formula is within 0.2% of our numerical estimates. Since the Hartree-Fock calculation 
is numerically limited to rather small skyrmion sizes,^! we believe that our field theoretical 
calculation incorporating short- distance physics will be a valuable tool to get a quantitative 
information over a wide range of skyrmion sizes. 



IV. SUMMARY 

We have obtained an effective Skyrmion action which incorporates short-distance physics 
as well as large-distance physics based on the Chern-Simons theory. We have numerically 
calculated the classical skyrmion profile and shown that for K > 2, our field theoretical re- 
sults exhibit an excellent agreement with the microscopic study. We have also demonstrated 
that the skyrmion energy very closely follows the formula suggested earlier by Sondhi et al. 
over a broad parameter range of interest. We believe that our field theoretical calculation 
incorporating short- distance physics will be a valuable tool to get a quantitative information 
for skyrmion over a wide parameter range of skyrmion sizes with reasonable numerical effort. 
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FIGURES 

FIG. 1. The number of spin flip as a function of g: The sohd curve corresponds to our field 

theoretical result and the open circles the Hartree-Fock calculation by Fertig et a/. [20]. 

FIG. 2. A(i^) as a function of the number of spin-flip K. The solid curve is calculated from 
the field theory. The filled circle represents the Hartre-Fock data and the open circles the exact 
diagonalization result by Abolfath et a/. [19]. 

FIG. 3. The skyrmion energy as a function of [g ln(^)]^/^. The line is a linear fit to our numerical 
data for small g. 
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